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Prufungsantrag gem. S 44 PatG ist gestellt 
© Verfahren und Einrichtung zur Raurnfilterung 

<§) Verfahren zur RaunrYfiftenjng elner Punkteverteilurtg, die 
cine von Stdrungen uberiagerte Struktur darstem und deren 
Punkte JeweUa durch n Koordi na ten werte einea n-dfmenaio- 
nalsn flaumes defJniert sind, baa welchem fur Jeden Punkt 
die Funktion, die die Abhangigkert der Anzahl der Punkte in 
elner Umgebung dea betreffenden Pun Jet as vom Durch mos- 
ser der Umgebung darateSt, sowie ein Skaltenmgskoeffi- 
adent, der gleioh dam Exponenten (a) einer eorfaehen 
Potenzfunktlon ist, weteha die genannte Funktion mdgOcriat 
gut annahert, ermittett warden, ferner die Skalierungskoeffi- 
zJerrtsn der Punkte elner der 2U fOtarnden Punktevertellung 
Squivalenten n-dimenalonaien stochaatischan Punfctevartej- ' 
lung ermJtteJt warden, die Olfferwu zwischen den Skulte- 
rungaJcoefnzientBn der zu ftitemden Punkteverteilung und 
m den SkaUerungskoeffiziefitan der st o ch aa ttechen Punkxever- 
tetiung ermittalt wird, und die Punkte, deren SkaJterungsko- 
J effizfenten in der Drfferenz enthattan dnd und die die 
weftgehend von Stoning an befreHe Struktur repreaentisren, 
in ainem mfndeatana zwekJtmenaionaJen Raum dargeatelh 
Oder weiterverarbeitet warden. 
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Beschreibung 

Die vorliegende Erfindung betrifft ein Verfahren und 
eine Binrichtung zur Raumfilterung, die es ermogiicheru 
Strukturen in n -dimensional en Raumen zu erkennen, 5 
attch wenn diese Strukturen stark durch Storungen 
Oberiagert sind 

Es ist bekannt, fur diesen Zweck lokale Sucfa verfah- 
ren for Dichtescttwankiingen, die sogenannte "Maxi- 
mum-Entrop ie- M ethode 1 * sowie die sogenannte °Maxi- 10 
mum-I jketihood-Metnode" zu verwenden. Diese be- 
kannten Verfahren erfordern Jedoch Annahmen oder 
Vorabinfonnationen urn Strukturen identifizieren zu 
kdnnen. Weitere Nachteile der bekannten Verfahren 
sind, daB eine Erweherung auf mehr als zwei Dimensio- 15 
nen sehr rechenaufwendig ist, daB uneinhehKche Mefi- 
gr6Ben und vieie un terschi edliche Korrelationen nicht 
verarbeitet werden kdnnen und daB die Empfindhchkeh 
bei irreguiftren Mustem mh starken unkorrelierten Sto- 
rungen gering ist 20 

Der voriiegenden Erfindung liegt dementsprechend 
die Aufgabe zugrunde, ein Verfahren und eine Einrich- 
tung anzugeben, mit denen auf verhfiltnismaBig einfache 
Weise n-dimenstonale Strukturen identifiziert werden 
kdnnen, die von starken Storungen Qberdeckt sind. 25 

Diese Aufgabe wird durch die in den Anspruchen 
gekennze ichnete und im folgenden nfiher erlauterte Br* 
findung gelosL 

GemfiB einem ersten Aspekt der Erfindung ist ein 
Verfahren zur Raumfilterung einer Ptmkteverteftung, 30 
die eine von Storungen tlberlagerte Struktur darstellt 
und deren Punkte jeweils durch n Koordinatenwerte 
eines n-dimensionalen Raumes definiert sind, gekenn- 
zeichnet durch die Verfahrensschritte: 

35 

a) fur jeden Punkt (Qm) wird die Funktion ermittelt, 
die die Abhangigkeit der Anzahl (N) der Punkte in 
einer Umgebung des betreffenden Punktes von ei- 
ner Abmessung (d) der Umgebung darstellt; 

b) fur jeden Punkt wird ein Skaherungskoeffizient 40 
(a(Qm)) ermittelt, der gleich dem Exponenten (a) 
einer einfachen Potenzfunktion 

N~d a 

ist, welche die im Schritt a) ermitteJte Funktion in 
einem gewissen Bereich der Abmessung (d) mSg- 
Iichst gut ann&hert, 

c) es werden die SkaUerungskoeffizienten der 
Punkte einer n-dimensionalen stochastischen Punk* 50 
teverteilung ernes der zu filternden Punktevertei- 
lung ftquivaienten Typs ermittelt, 

d) es wird die Differenz zwischen den Skalierungs- 
koefrizienten der zu filternden Punkteverteihmg 
und den SkaHerungskoeffizienten der stochasti- 55 
schen Punkteverteilung ermittelt, 

e) die Punkte, deren SkaHerungskoeffizienten in der 
Differenz enthaken sind und die die weitgehend 
von Storungen befreite Struktur reprasenueren. 
werden in einem cnindestens zweidimensionalen m 
Raum dargestelit 

GemaB einem zweiten Aspekt der Erfindung ist ein 
Verfahren zur Raumfilterung einer Punkteverteilung, 
die eine von Storungen Oberlagerte Struktur darstellt 65 
und deren Punkte jeweils durch n Kxrordinatenwerte 
ernes n-dimensionalen Raumes definiert sind, gekenn- 
zeichnet durch die Verfahrensschrhte; 
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a) fQr jeden Punkt (Qm) wird die Funktion ermittelt, 
die die Abhangigkeit der Anzahl (N) der Punkte in 
einer Umgebung des betreffenden Punktes von ei- 
ner Abmessung (d) der Umgebung darstellt; 

b) fur jeden Punkt wird ein Skalierungskoeffizient 
(a(Qm)) ermittelt, der gleich dem Exponenten (a) 
einer einfachen Potenzfunktion 

N~d» 

ist, welche die im Schritt a) ermitteJte Funktion in 
einem gewissen Bereich der Abmessung (d) mog- 
tichst gut ann&hert, 

c) es werden die Skafierungskoeffizienten der 
Punkte einer n-dimensionalen stochastischen Punk- 
teverteilung eines der zu filternden Punktevertei- 
lung Squivalenten Typs ermittelt, 

d) jedem vorkommenden Wert der SkaHerungsko- 
effizienten der zu filternden und der stochastischen 
Punkteverteilung wird eine Kumulativzahl (P(< a)) 
zugeordnet, die gleich der Anzahl der ermittelten 
SkaherungskoefTrzienten der betreffenden Vertei- 
iung ist, die grfiBer oder gleich dem betreffenden 
SkaUerungskoeffizienten sind, dividiert durch die 
Cesamtzahl der gemessenen Skalierungskoern- 
zienten der betreffenden Verteiiung, wobei die in 
diesem Schritt ermittelten Kumulativzahl en der 
stochas t ischen Verteiiung so normiert werden, daB 
sie bei den grdBten Zahlenwerten mit den grofiten 
Zahlenwerten fur die zu ftltemde Verteiiung Ober- 
emstunmen, 

e) die Wertebereicfae der Skaliexungskoeffizienten 
der zu filternden und der stochastischen Punktever- 
teilung werden in TeObereiche aufgeteilt; 

f) fur jeden Tedi bereich wird die Differenz der zuge- 
horigen. im Schritt d) ermittelten Kumulativzahlen 
der zu filternden Punkteverteilung und der stocha- 
stischen Punkteverteilung erzeugt, und 

g) die Punkte, deren SkaHerungskoeffizienten in 
der Differenz enthalten sind und die die weitgehend 
von Storungen befreite Struktur reprasentieren, 
werden in einem mindestens zweidimensionalen 
Raum dargestelit 

Vorzugsweise umfafit der Schritt a) folgende Teil- 
schritte: 

al) die Punkteverteilung wird abgetastet und die 
da bei ermittelten Koordinatenwerte der einzelnen 
Punkte (Qm) werden gespekhert; 
a2) fur jeden Punkt werden Koordinatenwerte, die 
den Koordinaten des Punktes benachbart sind, in 
sukzessive grdBeren Abstfinden durch Durchsu- 
chen des Speichers abgefragt, ob diese Koordina- 
ten einen Punkt enthahen und die Anzahl der dabei 
ermittelten Punkte wird als Funktion des Abstan- 
des gespeichert 

Der Verf ahrensscfaritt b) umfaBt ferner vorzugsweise 
die folgenden Teilschritte: 

bl) es wird die Funktion log N~a log d ermittelt, 
wobei N die Anzahl der Punkte und d eine Dimen- 
sion des den betrachteten Punkt umgebenden Be- 
reiches darstellt, und 

b2) es wird die Steigung eines im wesentlichen ge- 
raden Tefles dieser Funktion ermittelt und als Ska- 
lierungskoefflzient verwendet. 
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Erne bevorzugte Einrichtung zur Durchf Qhrung eines 
Verf aureus gem&B der Erfindung enth&it: 

a) eine Einrichtung (100) zum Ermittein der Koordi- 
natenwerte jedes Punktes; 5 

b) einen Speicher (102) zum Speichern do* Koordi- 
natenwerte jedes Punktes; 

c) eine Rechenemricbtung (104, 10& 108) die fQr 
jeden Punkt einen Skafierungskoeffizieiiten ermit- 
telt, der gleich dem Expo aen ten (a) einer einfachen io 
Poteiizfunktion 

N~d* 

ist, wekhe die Abhfingigkeit der Anzahl (N) der 15 
Punkte in einer Umgebung des betreffenden Punk- 
tes von einer Abmessung (d) der Umgebung dar- 
steiit; 

d) einer Vorricfatung (112) zum BereitsteDen von 
Skaflerungskoeffizienten der Punkte einer n-di- 20 
mensionaien stochastischen Punkteverteilung eines 
der zu filternden Punkteverteilung aqirivalenten 
Typs; 

e) eine Vorrichtung (116) zum Selektieren der 
Punkte der zu filternden Punkteverteilung, die in 23 
exnem Bereich der Skalierungskoeffizienten liegen, 

in denen sich die Skalierungskoeffizienten der zu 
filternden Punkteverteilung wesentllch von den 
Skalierungskoeffizienten der stochastischen Punk- 
teverteilung unterscheiden.und 30 

f) eine Vorrichtung zur Darstelhing der Punkte, de- 
nen die selektierten Skalierungskoeffizienten zuge- 
ordnet sind 

Das Verfahren und die Einrichtung geraaB der Erfin- 35 
dung sind dem Stand der Technik erheblkh uberiegen 
und erlauben die Identifizierung von Strukturen, die 
durch so starke Storungen verdeckt sind, dafi die be- 
kannten Verfahren versagen. Der erforderliche zeitiiche . 
und apparative Aufwand ist gering. 40 

Das Verfahren und die Einrichtung gemaB der Erfin- 
dung sind auf Wertetupei beliebiger Dimensionen an- 
wendbar. Bevorzugte Anwendungsgebiete sind die au~ 
tomatische Bfldauswertung, z. B. in der Medizin, bei Sa- 
tellhenbildern und Luftauf nahmen sowie bei der Werk- 45 
stoffprufung. 

Das Verfahren und die Einrichtung gemaB der Erfin- 
dung lassen sich auf beiiebige Punktmengen in n-dimen- 
sionalen Kooniinatenraumen (n — naturliche Zahl) an- 
wenden. Dabei kann der BegrifT 'Punkf sowohl im ma- so 
thematischen Shone verstanden werden ais auch als Ko- 
ordinatentupel zur IdenuEkation von Bildpixeln und 
dergieichen. Die entsprechenden Koordlnatenrflume 
kannen beiiebige gauzzahlige Dimensionen haben, die 
in der Praxis mindestens gleich 2 sein werden. Die Koor- 55 
dinaten sind nicht auf eine einheitltche physikahsche Di- 
mension beschrankt, sondern konnen z. B. Ortskoordi- 
naten sowie einen oder mehrere zusatzliche Parameter, 
wie einen Farbwert, einen Energiewert oder dergiei- 
chen umfassen. In einem drekiimensionaleu Koordina- eo 
tenraum k6nnen beispieisweise die ersten betden Koor- 
dinaten den Ort in einer Flache und die dritte Koordina- 
te einen auf eine geeignete Skala transformierten Ener- 
gies oder Farbwert enthahen. Im einfachsten Falle wtrd 
man die Skalierung der Koordtnaten so wahien, dafl die 65 
resulderende Punktmenge in ein Einhextsquadrat oder 
einen Einhehswurf el abgebildet werden kann. 

Bei einer sich zehlich todernden MefigrOBe fZeitrei- 



he*) wird ein n-dimensionaler Koordinatenraum ge- 
wthlt, wobei die Punktkoordinaten durch n aufeman- 
derf olgende Messungen der interess i erenden Gr&Be be- 
snxnmt sind. 

Im folgenden wird die Erfindung unter Bezugnahme 
auf die Zeichnungen naher eriaiitert, dabei werden noch 
weitere Vorteiie und Merkmaie der Erfindung zur Spra- 
chekommen. 

Eszeigen: 

Fig. la eine veremfachte Daxsteflung einer in erster 
Naherung einem PunM entsprechenden (^0-dimensiona- 
len") Verteilung von schwarzen Pixeln; 

Fig* lb ein Diagramm, das schema tisch die Anzahl N 
der schwarzen Pixel inner halb eines als krasfdrmig an- 
genommenen Umgebungsgereiches urn ein vorgegebe- 
nes, durch einen ausgezogenen Kreis bezeichnetes Pixel 
in Abhangigkeh vom Radius des Umgebungsbereiches 
darsteflt; 

Fig. 2a und 2b Darsteflungen entsprechend Fig. la 
bzw. 1 b fur eine Iineare (l-dimensk>nak) Verteilung von 
schwarzen Pixehv, 

Fig. 3a und 3b Darsteflungen entsprechend Fig. la 
bzw, lb fur eine flfichige (2-dimensronale) Verteilung 
von schwarzen Pixeln; 

Fig. 4a eine zweidimensionale Punktmenge mit 900 
stochastisch verteilten schwarzen Punkten auf weiBem 
Untergrund; 

Fig. 4b die Punktmenge gema\B Fig. 4a in deren Mitte 
zusatzlich 100 Punkte angeordnet sind, die einer Struk- 
tur mit GauBverteilung entspricht; 

Fig. 4c das Spektrum von Skalierungskoeffizienten 
der Punktverteflung gemfiB Kg. 4a; 

Fig, 4d das Spektrum der Skalierungskoeffizienten 
der Punktverteflung gemaB Fig, 4b; 

Fig. 4e die kumulative Verteilungsfunktion fur die 
Punktverteflung gemaB Fig- 4a; 

Fig. 4f die kumulative Verteilungsfunktion der Punkt- 
verteflung gemaB Fig. 4b; 

Fig. 5a eine nut dem Rdntgensatclliten ROSAT ge- 
messene Photonenkarte des Ousters A 1314; 

Fig. 5b das Spektrum der SkalierungskoefTodenten 
fur die Punktverteflung gemaB Fig. 5a; 

Fig. 5c die kumulative Verteiajngsfunktion fur die 
Punktverteflung gemflB Fig. 5a; 

Fig. 5d die Abweidumg der kumuiauven Verteilungs- 
funktion gemaB Fig. 5c von der einer statistischen 
Punktverteflung oder -menge; 

Fig. 6a ein Spektrum der Skalierungskoeffizienten f Or 
die Punkteverteilung gemaB Fig. 5a unter Hinzunahme 
einer Energiekoordinate; 

Fig. 6b eine Darsteflimg der Ortskoordinaten und der 
Energie der Punkte, die zu dem schraffierten Bereich 
des Spektrums der Skalierungskoeffizienten gemaB 
Fig. 6a nach Skalierung auf einem Emheitswurfel geh&- 
ren; 

Fig. 6c die Projektion der Verteilung im Einheit3w0r- 
fei auf die X- Y-Ebene; 

Fig. 7a eine Punktverteflung, welche einem von Rau- 
schen Qberiagerten Sinussignai entspricht; 

Fig. 7b ein Spektrum der Skalierungskoeffizienten 
der von der stochastischen Stdrung Qberiagerten Sinus- 
funktton fur einen fanfdimensioaalen Koordinaten- 
raum; 

Fig. 7c ein Spektrum der Skalierungkoeffizienten der 
von der stochastischen St&rung Qberiagerten Sinusfunk- 
don fur einen si ebendhn ensionalen Koordina tenraum ; 

Fig. 7d eine Darsteflung der rekonstruierten Struktur, 
die aus den Ennittlungsskalierungskoeffizienten der 
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vorgegebenen Verteilung abzuglicn der Skalierungsko- 
efflezieiiten fur eine stochastische Verteilung ermitteit 
wurde; 

Fig, 8 ein ScfaaJtbild einer Flnricfatung zur Durchfuh- 
rung des voriiegenden Verfahrens* 

lm f olgenden soil der Begriff *Punkt v sowofal im ma- 
thematischen Sinne verstanden werden als auch Koor- 
dinatentupel bedeuten, die jewcils ein Bildpixel o. dgL 
definieren. Die Koordmatenrauine, in die diese Punkte 
Qm eingebettet sand, konnen beliebtge ganzzahlige Di- 
mensionen, im ailgemeinen 2: 2» annehmen. Die Koordi- 
naten sind. wie erwahnt tticfat auf erne einheltliche phy- 
rikalische Dimension, z, B. Lange, beschrankt 

Der erste Schritt des voiiiegenden Verfahrens be- 
steht darin, fur jeden Punkt der zu imtersuchenden 
Punkteverteilung oder Punktmenge (oder, was Fflr den 
zweidimensionalen Fall Equivalent ist fur jedes Pixel 
eines Rasters oder xy-Feldes) einen sog. Skalierungsko- 
effizienten zu bestunmen. Der Skaiierungskoeffizient 
wild aus der Punktdkhte in der Umgebimg des betrach- 
teten Punktes, d h. der Anzahl N der Punkte im jeweili- 
gen Umgebungsbereich, ermittelL Insbesoudere wird 
unter der Annahme, daB die Punkteanzahl N in der 
Umgebung des betrachteten Punktes in Abhangigkert 
vom Durchmesser d dieser Umgebung einer einfachen 
(d. h. nur ein emziges Potenzglied enthaltenden) Potenz- 
funktion 

N = d a 

folgt der Skaiierungskoeffizient fur diesen Punkt als der 
Exponent a der Potenzf unktion defimert 

Ob die Punktanzahl In einer beschrankten Umgebung 
eines betrachteten Punktes in etwa einer einfachen Po- 
tenzfunktion folgt, und der Exponent dieser Potenz- 
f unktion, lassen sich am einfachsten aus einer doppeltlo- 
garithmiscben Darstellung der Punkteanzahl N als 
Funktion des Durchmessers d ermitteln: 

logN-alogcL 

Bei einer solchen Darstellung ergibt sich im Falle ei- 
ner einfachen Potenzf unktion eine Gerade. Die Stei- 
gung dieser Geraden ist dann der gesuchte Skaiierungs- 
koeffizient fur den betrachteten Punkt Bei nichtkreis- 
f Ormiger Umgebung, z. B» im Falle einer quadra tischen 
Umgebung eines Pixels eines xy-Rasters, bedeutet d ei- 
ne andere entsprecheode Abmessungsdimension, z. B. 
die Seitenlange ernes die Umgebung biidenden Qua- 
drats, 

Fur jeden Punkt Q der Punktmenge werden der zuge- 
hdrige Skaiierungskoeffizient bq und die zugehdrigen 
Punktkoordinaten xq, y<j, zq, • . . gespeichert 

Nachdem die Skahenmgskoeffizienten fur alle Punk- 
te der zu untersuchenden Punktmenge ermitteit warden 
sind, wird die rUufigkehsverteilung der Skalierungsko- 
effizienten ermitteit, d. fa- die Anzahl f(a) der jeweiligen 
Skalierungskoeffizienten (a) in Abhftngigkeit von ihrem 
Wert Zusatzlicb wird gespeichert welche Punkte der 
untersuchten Punktmenge zu den jeweiligen Werten 
der Haufigkeitsvertetlung der Skalierungskoeffizienten 
gehdren. 

Fur eine vollkommen statistische Punktmenge im 
Grenzfall zum Koutinuum haben alle Punkte den glei- 
chen Skalierungskoeffizienten. Die Haufigkehsvertei- 
lung reduziert sich dann also auf einen einzigen Wert 
Diese Haungkehsverteilung wird nun so normiert, daB 
in einem n-dimensionalen Koordinatenraum die Skalie- 
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rungskoeffizienten alle den Wert n annehmen und der 
Wert der Hfiufigkei tsverteilung ebenfalls n betrfigt 

Die transformierte Haufigkeitsverteflung wird als 
Spektrum der SkaHenmgskoefixzdenten bezeichnet Das 
5 Spektrum der Skalierungskoeffizienten ist eine eindeu- 
tige Abbfldung der lokalen und globalen Egenschaften 
der untersuchten Punktmenge. 

Analog zu dem obigen Vorgehen wird ferner das 
Spektrum der Skaiierungsko effbden ten einer der unter- 
io suchten Punktmenge entsprechenden Punktmenge (also 
z. B. einer unter den gtekhen oder analogen Bedingun- 
gen gewonnene Punktmenge) bestimmt die keine Dicfa- 
te- oder Korrelationsgradienten aufweist Dieses Spek- 
trum kann entweder auf analytischem Wege gewonnen 
15 werden, wie bei dem unten angefuhrten Beispiel, oder 
durch Simulation bzw. geeignete Messung einer ent- 
sprechenden Punktmenge. Fur diese stochastische Ver- 
gieichsmenge ist es nicht notwendig, die Punkte den 
jeweiligen Werten der HAufigkehsverteilung zuzuord- 
20 nen, es wird nur das Spektrum der Skalierungskoeffi- 
zienten als sokhes benotigt 

Es wird nun der Unterschied oder die Abweichung 
des SkaHenmgskoen^enten-Spektrums der untersuch- 
ten Punktmenge vom SkaHerungskoeffizienten-Spek- 
2s trum der Vergleichsmenge ermitteit Der Grad der Ab- 
weichung dieser Spektren, der den dichten bzw. stark 
korrelierten Gebieten entspricht ist ein Mafi fur das 
Auftreten von Dichte bzw. Korrelationsgradienten* Die 
Punkte der untersuchten Menge, die diesem abweichen- 
30 den Tefl zugeordnet sind, identifizieren die °Struktur w 1 
die die Abweichung der Dichte bzw. Korrelationsgra- 
dienten bedingt 

Die Abweichung kann auf zwei Art en bestimmt wer- 
denz Bei der ersten Methode werden, wie oben beschrie- 
35 ben, die Spektren der Skalierungskoeffizienten mhen> 
ander verglichen. Urn zur Haungkehsverteilung zu ge- 
langen ist es in praktischen Anwendungen notwendig, 
die gefundenen Skalierungskoeffizienten in Wertekana- 
le oder Wertebereiche, sog. "Bins", einzuteOen. Durch 
40 die Breite der Wertekanfile wird der statistischc Fehler 
fur die Werte der HSufigkeitsverteilung bestimmt Eine 
groBe Wertekanalbrerte bedeutet einen kleinen stausti- 
schen Fehler, allerdings auch eine geringe Aufldsung. 
Bei der zweiten Methode zur Besthnmung der Ab- 
45 weichung von gegebener zu erwarteter Verteilung wird 
die EmfQhnmg eines zusatzlich en Parameters in Form 
der Binbreite vermieden. Statt des Spektrums der Ska- 
lierungskoeffizienten wird die kumulative oder Integra- 
le Verteihmgsfunktion der Skalierungskoeffizienten be- 
so stimmt Dazu wird jedem vorkommenden Wert der Ska- 
lierungskoeffizienten die Zahl der gemessenen Skalie- 
rungskoeffizienten zugeordnet, die grd&er oder gleich 
diesem Wert sind. Diese Zahl wird durch die Gesamt- 
zahl der gemessenen Skalierungskoeffizienten dividiert 
55 Die Abhangigkeit der bei der Division erhaltenen Werte 
(P(a)) vom Wert (a) des Skalierungskoeffizienten stelit 
die kumulative Verteilungsfunktion dar. Analog wird die 
kumulative Verteilungsfunktion fur eine Vergleichs- 
menge bestimmt. Die kumulative Verteilungsfunktion 
eo der Vergleichsmenge wird so normiert daB sie rait der 
kumulativen Verteilungsfunktion der untersuchten 
Menge bei den groBten vorkommenden Werten der 
Skalierungskoeffizienten flbercinstimmt Analog zu den 
Spektren ist die Abweichung der kumulativen Vertei- 
65 lungsfunktionen bei kleinen Werten der Skalierungsko- 
effizienten ein MaB fflr das Auftreten von Dichte- und 
Korrelationsgradienten, 

Bevor auf spezielle AusfQhrungsbeispiele eingegan- 
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gen wird, soil das der Erfindung ztignmdeliegende Prin- 
jxp unter Bezugnahme auf die Fig. 1 bis 3 erlautert wer- 
den.Dk Fig. la, 2a und 3a stelkn jeweils ein zweidimen- 
sionales Pixel-Raster dar. "Scfawarze" Pixel, die •Punkte^ 
auf einem weiBen Hintergrund darsteilen, sind durch ein 
Kreuz bezeichnet 

Fig. 1 zeigt eine fleckartige, in erster N aliening 
p unktar tige oder nuHdxmensionale Verteihing von 
schwarzen Pixeln im Zentnxm des Rasters. Die Ermitt- 
lung des Skalienmgskoeffizienten sei am Beispiel des 
mittieren sciiwarzen Pixels, das durch einen ausgezoge- 
nen Kreis bezeidinet ist eriauiert Fur dieses Pixel wird 
die "Dichte", genauer gesagt die Anzabl der Punkte in 
Abhangigkeit vom Durchmesser d des das Pixel umge- 
benden Bereichs (oder vom Abstand r von der Mitte des 
betreffenden Pixels aus gerecbnet) bestnnmt Wie in 
Fig. lb dargestellt ist, steigt die Anzahi der schwarzen 
Pixel im betrachteten Bereich, mit zunebmendem d zu- 
erst an, bleibt dann aber konstant Der Grad der Potenz- 
funktion, die den konstanten Hauptteil der Kurve in 
Fig; lb darsteilt, und damit der Skalierungskoefnzient 
(oder Skalierungsexponent) der Dicbte-Abstand-Funk- 
tion ist dementsprechend null 

Fig. 2a zeigt eine lineare Verteihing von schwarzen 
Pixeln. Auch hier sei das durch einen ausgezogenen 
Kreis bezeicbnete Pixel betrachtet und es wird die An- 
zahi der schwarzen Pixel in der Umgebung dieses Pixels 
ais Funktion des Durchmessers bestimmt wie durch ge- 
stricheite Kreise angedeutet ist Hier nimmt offensieht- 
lich die Anzahi der Pixel linear mh dem Durchmesser d 
zu. Die Dichte-Durchmesser-Funktion ist eine Gerade 
und der Skalieningskoeffizient ist dementsprechend 
gietcheins. 

Fig. 3a zeigt eine fUchige Verteiiung von schwarzen 
Pixeln. Hier nimmt die Anzahi der schwarzen Pixel pro- 
portional zum Quadrat des Durchraessers d zu, Der 
Skalieningskoeffizient ist dementsprechend zweL 

Fur Verteilungen, die von den oben bcschriebencn 
abweichen, ergeben sich andere, ggf. nicfat-ganzzahlige 
Skaherungskoeffeienten. 

Fig. 4 zeigt die hohe Nachweisempfhidlichkeit des 
vorliegenden Verfahrens am Beispiel einer Ponktcnenge 
in einem zweidimensionalen Koordinatenraum. In 
Fig. 4a ist eine stochastische Punktmenge mit 900 Punk- 
ten dargestellt Fig. 4b zeigt die Punktmenge gemaB 
Fig. 4a, die in der Mitte zusatzlich 100 Punkte mit einer 
gauflformigen Dichteverteilung (°Struktur s ) enthalt 
Der la-Radius dieser GauBverteilung betrdgt 1/8 des 
Felddurchmessers. Die Struktur befindet sich innerhalb 
des eingezeichneten Kreises. 

Fig. 4c und Fig. 4d zeigen das Spektrum der Skalie- 
rungskoeffizienten fur die stochastische Punktmenge 
gemaB Fig. 4a bzw. for die die Struktur enthaltende 
Punktmenge gemafl Fig. 4b. In den Fig. 4e und 4f ist die 
kumulative Verteflungsfunktion fflr die Punktmenge ge- 
maB Fig. 4a bzw. 4b dargestellt, 

Bei der PrOfung, ob eine gegebene Punktmenge eine 
Struktur enthalt wird man im allgemeinen die kum ulati- 
ve oder integrate VerteOungsfunktion der zu untetsu- 
chenden Punktmenge an die theoretisch erwartete Kur- 
ve fur eine entsprechende stochastische Punktmenge 
anpassen und aus dieser Anpassung die Anzahi der im 
untersuchten Feld stochastisch verteilten Punkte be- 
sthnmen. Weicht die so bestimmte Zahl stochastiscner 
Punkte signifikant von der Zahl der im untersuchten 
Feld insgesamt vorhandenen Punkte ab, so enthalt das 
untersuchte Feld eine Struktur, wie sie oben definiert 
wurde. Wendet man dieses Krtterium auf 100 statistisch 



unabhSngige Felder an, die gemaB dem im vorherge- 
henden Abschnht beschxiebenen Verfahren erzeugt 
wurden, so werden bei einer Abweichung von mehr ais 2 
a 95 Felder ais nichtstochastisch identifiziert wobei in 
5 88 Fallen die Zahl der stoachastischen Punkte (900 
Punkte) bis auf la (30 Punkte) genau bestimmt wird. 

Die kumulative VerteOungsfunktion laSt also erken- 
nen, ob die betrachtete PunkteverteOung eine Struktur 
enthalt 

1Q Ais prakiisches Beispiel zeigt Fig. 5a eine mit dem 
Rdntgensatelliten ROSAT gemessene Photonenkarte 
des Clusters A1314. Jedem vom Detektor des Satelfiten 
registrierten Photon sind zwei Ortskoordinatenwerte 
und ein Energiewert zugeordnet In Fig- 5a entspricht 
15 jeder schwarze Punkt in dem zweidimensionalen Koor- 
dinatenraum einem gemessenen Photon, das durch seine 
Detektorkoordinaten identifiziert wird. Fur diese 
Punktmenge wind das Spektrum der Skalierungskoeffi- 
zienten bestimmt das in Fig. 5b dargestellt ist Dieses 
20 Spektrum wird mit dem Spektrum einer Menge vergti- 
cben, deren Punkte stochastisch verteilt sind. FOr eine 
solche zufalfig ausgespielte, stochastische Punktvertei- 
lung lafit sich das Spektrum der Skafierungskoeffizien- 
ten analytisch bestimmen. 
25 Fig. 5c zeigt die kumulative Verteilungsfunktion der 
Punktmenge gemaB Fig. 5a (Kreuze) und die kumulati- 
ve Verteilungsfunktion der stochastischen Punktmenge 
(ausgezogene Kurve). 

In Fig. 5d ist die Abweichung der kumuktiven Vertei- 
30 lungsfunktion der Punktmenge gemaB Fig. 5a von der 
kumulativen Verteilungsfunktion der stochastischen 
Menge dargestellt Die Abweichung wurde im vorlie- 
genden Fafie noch dadurch verstarkt dafl man die Ener- 
gie ais dritte Koordinate hinzunahm. Die Energie und 
35 die Ortskoordinaten wurden auf einen Einheitswurfel 
skaliert und von der so erhaltenen Punktmenge das 
Spektrum der SkaHerungskoeffizienten bestimmt das in 
Fig. 6a dargestellt ist Der Teil des Skalierungskoeffi- 
zienten-Spektrums, der auf die die Abweichung verur- 
40 sachenden Punkte in Fig. 5a zurflckgeht ist in Fig. 6a 
schraffiert dargestellt Stellt man diese Punkte in einem 
Einheitswurfel dar, so erhah man die in Fig. 6b darge- 
stellte raumhche Punktevertetlung. 

Fig. 6c zeigt die Projektion dieser Punktverteilung in 
45 die xy-Ebene in der auch die ursprQngUche Punktever- 
teilung gemaB Fig. 5a dargestellt ist Die durch das vor- 
hegende Verfahren identifizierte Struktur gemaB 
Fig. 6c gibt die Position und Konfiguration des Clusters 
A1314 sehr gut wieder. 
so Ais weiteres Beispiel zeigen die Fig. 7a bis 7d die 
Auswertung einer PunkteverteQung (Fig. 7a) entspre- 
chend einem von Rauschen uberiagertem unvoflstandi- 
gen SinussignaL Aus einer shnuherten Zeitreihe, die sich 
aus einer unvollstandigen Sinuskomponente und einer 
55 vollstfindig stochastischen Komponente zusammensetz- 
te, wurden n-dimenaionale Koordinatenraiime konstru- 
iert (n = 2, 5, 7> Jeder Punkt in Fig. 7a ist durch n 
aufeinanderfolgende MeBpunkte definiert, wobei die er- 
ste Koordinate der Punkte die gesamte Zeitreihe durch- 
60 lauft Ais Vergieichsmenge wird eine Zeitreihe verwen- 
det die nur eine stochastische Komponente enthalt Wie 
, die Fig. 7b und 7c far n « Sbzw.n = 7 zeigen, wird die 
Trennung der stochastischen Komponente von der si- 
nusartigen korrelierten Komponente mit zunehmender 
65 Dimension des Koordinatenraumes groSer. 

Fig. 7d zeigt cue Rekonstruktion des Sinussignals aus 
einem Skah'erungskoeffizienten-Spektrum fur n =» 2. ^ 
Fig. 8 zeigt ein Funktions-Blockschaltbild einer Ein- 



DE 43 17 

9 

richturtg zur Durchfahrung des oben beschriebenen 
RaumfUterungsverfahrens auf der Basis der Skalie- 
niiigskoeffizientensp^ktren. Die einzelnen Emheiten 
der Einrichtung sind dabei durch ihre Funktionen be- 
zeichnet 5 

Die zu ftkemde Punkteverteilung wind als ernes in 
einer Abtaatemrichtung 100 zum Ermitteln der Koordi- 
nfitenwerte der Punkte Qm der Punkteverteilung abge- 
tastet. 

Die ermittelten Koordinatenwerte XQm, yQm, . . . wer- io 
den in einer Spekhereinheit 102 gespeichert 

An die Speichereinheit 102 isl eine Integriereinheit 
104 angescbiossen, welcfae die in der Speichereinheit 
102 gespeicherten Punktekoordmaten abfragt und die 
Punkteanzahi/Bereicfasgr6Be-Funktion N (d) ermittelt is 

Eine Logarithzniereinheit 106 ermittelt aus der 
Punkteanzafal/BereichsgroBe-Funktion des betreffen- 
den Punktes den SkaUerungskoeffizienten a(Qm) des be- ' 
treff enden Punktes Qm auf der Basis der f olgenden For- 
mel 20 

log (N2 - Nl) 
log (£2 - di) 

wobei Nl und N2 die Aozahlen der Punkte in den Beret- 
chen mit den Abmessungen dl bzw. 62 bedeuten. Man 
erhait dadurch den Skaiierungskoeffizienten a (Qm) des 
betreffenden Punktes der mh den Punktekoordinaten in 30 
einem Speicher 108 gespeichert wird Man kann hierzu 
den Speicher 102 verwenden und den schon gespeicher- 
ten Koordinatenwerten des betreffenden Punktes den 
ermittelten Skalierungskoeffizienten zuordnen. 

Als n&chstes werden die Skalierungskoeffizienten aus 3s 
dem Speicher 108 abgefragt und in einer Diskrinunato- 
reinheit 110 in Wertebereiche eingeordneL 

Die Einrichtung gem. Fig. 8 enthSlt ferner einen Spei- 
cher 112, der dem Speicher 102 entspricht und die Koor- 
dinatenwerte einer mit der auszuwertenden Punktever- 40 
teilung vergleichbaren stochastischen Punkteverteilung 
enth&lt. In einer an dieseu Speicher angeschiossenen 
Einheit 1 14 werden die Skalierungskoeffizienten fur die- 
se stochastische Verteilung ermittelt und in Werteberei- 
che entsprechend den Werteberekhen der Einheit 110 45 
eingeordneL 

In einer Differenzstufe 116 werden die Skalierungs- 
koeffizienten der stochastischen Verteilung von den 
Skalierungskoeffizienten der auszuwertenden Punkte- 
verteilung abgezogen. SchlieBlich werden die Punkte, so 
die den bei der Differenzbildung verbliebenen Skalie- 
rungskoeffizienten entsprechen, in einer Ausgabeein- 
heit 1 18 *ruckprojiziert B , & h. sie werden in einem n-di- 
mensionalen Raum wie die abgetastete Punktevertei- 
lung mittels ernes Druckers und/oder eines Monitors 55 
dargesteUt oder als gefiherte Werteverteilung in ande- 
rerWeise weiterverarbeitet. 

Die Stufen 102 bis 116 kdnnen durch diskrete Bauein- 
hehen oder durch einen Mikroprozessor realisiert wer- 
den* 60 

Patentanspriiche 

1. Verf ahren zur Raumfilterung einer Punktevertei- 
lung, die eine von Stdrungen uberlagerte Struktur 65 
darsteQt und deren Punkte jewefls durch n Koordi- 
natenwerte eines n-dimensionalen Raumes defi- 
niert And. gekennzeichnet durch die Verfahrens- 
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schritte: 

a) fur jeden Punkt (Qm) wird die Funktion er- 
mittelt, die die Abbfingigkeit der Anzahl (N) 
der Punkte in einer Umgebung des betreffen- 
den Punktes von einer Abmessung (d) der Um- 
gebung darsteilt; 

b) fur jeden Punkt wird ein Skaherungskoeffi- 
zient (a(Qm)) ermittelt. der gleich dem Expo- 
nenten (a) einer einfachen Potenzfunktion 

N~d a 

ist, welche die im Schritt a) ermittelte Funktion 
in einem gewissen Bereich der Abmessung (d) 

mflgBcfaSt g»t annBher t, 

c) es werden die Skalierungskoeffizienten der 
Punkte einer n-dimensionalen stochastischen 
Punkteverteilung eines der zu fUternden Punk- 
teverteilung Iquivalenten Typs ermittelt, 

d) es wird die Differenz zwischen den Skalie- 
rungskoeffizienten der zu fihernden Punkte- 
verteilung und den Skalierungskoeffizienten 
der stochastischen Punkteverteilung ermittelt, 

e) die Punkte, deren Skalierungskoeffizienten 
in der Differenz enthalten sind und die die 
weitgehend von Stdrungen befreite Struktur 
reprSsentieren, werden in einem mindestens 
zweidimensbnalen Raum dargestellL 

2. Verf ahren zur Raumfilterung einer Punktevertei- 
lung, die eine von Stdrungen flberiagerte Struktur 
darsteilt und deren Punkte jeweils durch n Koordi- 
natenwerte eines n-dimensionalen Raumes defi~ 
niert sind, gekennzeichnet durch die Verfahrens- 
schritte: 

a) fQr jeden Punkt (Qm) wird die Funktion er- 
mittelt, die die Abhlngigkett der Anzahl (N) 
der Punkte in einer Umgebung des betreffen- 
den Punktes von einer Abmessung (d) der Um- 
gebung darsteilt; 

b) f Or jeden Punkt wird ein Skafierungskoeffi- 
zient (a(Qm)) ermittelt, der gleich dem Expo- 
nents (a) einer einfachen Potenzfunktion 

N~d* 

ist, welche die im Schritt a) ermittelte Funktion 
in einem gewissen Bereich der Abmessung (d) 
moglichst gut ann&hert, 

c) es werden die Skalierungskoeffizienten der 
Punkte einer n-dimensionalen stochastischen 
Punkteverteilung eines der zu fihernden Punk- 
teverteilung iquivalenten Typs ermittelt, 

d) jedem vorkommenden Wert der Skalie- 
rungskoeffizienten der zu fUternden und der 
stochastischen Punkteverteilung wird eine Ku- 
mulativzahl (P(a)) zugeordnet, die gleich der 
Anzahl der ermittelten Skalierungskoeffizien- 
ten der betreffenden Verteilung ist, die grdBer 
oder gleich dem betreffenden Skalierungsko- 
effizienten sind, divtdiert durch die Gesamt- 
zahl der gemessenen Skalierungskoeffizienten 
der betreffenden Verteilung, wobei die in die- 
sero Schritt ermittelten Kunniiatrvzahlen der 
stochastischen Verteilung so normiert werden, 
daS sie bei den groBten Zahlenwerten mh den 
gro&ten Zahlenwerten fur die zu filternde Ver- 
teilung Obereinstimmen, 

c) die Wertebereiche der Skalierungskoefri- 
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zienten der zu fiUemdec und dcr stochasti- 
schen Punkteverteilung werden in Teilberei- 
che auf geteilt; 

f) fflr jeden Teilbereicb wird die Differenz der 
zugehfcrigen, im Schritt d) ennitteiten Kumu- 5 
lativzahlen der zu filteroden Punkteverteilung 
and der stochastischen Punkteverteihmg er- 
zeugtund m 

g) die Punkte, deren Skaherungskoeffizienten 

in der Differenz enthahen mud and die die 10 
wehgehend von StOnmgen befreite Struktur 
reprfisentieren, werden in einem mindestens 
zweidimensionalen Raum dargesteUt 

3. Verfahren nach Anspruch 1 oder 2, dadurch ge- 
kennzeichnet, daB der Scfarht a) folgende Teil- 15 
schritte umfafit: 

al) die Punkteverteilung wird abgetastet und 
die dabei ennitteiten Koordmatenwerte der 
einzelnen Punkte (Qm) werden gespeichert; 
a2) fUr jeden Punkt werden Koordinaten werte, 20 
die den koordinaten des Punktes benachbart 
sind, in sukzessive gr&Beren AbstSnden durch 
Durchsuchen des Speichers abgefragt, ob die- 
se Koordinaten einen Punkt enthalten und die 
Anzahl der dabei ennitteiten Punkte wird als 25 
Funktion des Abstandes gespeichert 

4. Verfahren nach einem der Anspruche 1 bis 3, 
dadurch gekennzeichnet, daB der Verfahrensschritt 
b) die foigenden Teilschritte enthaJt: 

bl) es wird die Funktion log N — a log d ermh- 30 
telt, wobei N die Anzahl der Punkte und d eine 
Dimension des den betrachteten Punkt umge- 
bendcn Bereiches darstellt, und 
b2) es wird die Steigung eines im wesentiichen 
geraden Teiles dieser Funktion ermittelt und 35 
als Skalie rungs koeffizient verwendet- 

5. Hinrichtung zur Raumfilterung einer Punktever- 
teilung, die eine von Storungen Oberlagerte Struk- 
tur darstellt und deren Punkte jewefls durch N Ko- 
ordmatenwerte eines n-dimensionalen Raumes de- 40 
finiert sind, gekennzeichnet durch 

a) eine Einricbtung (100) zum Ermittein der 
Koordinatenwerte jedes Punktes; 

b) einen Speicher (102) zum Speichem der Ko- 
ordinatenwerte jedes Punktes; 45 

c) eine Recheneinricfatung (104, 106, 108) die 
fiir jeden Punkt einen Skalierungskoeffizien- 
ten ermittelt, der gleich dem Exponenten (a) 
einer einf achen Potenzfunktion 

50 

N~d a 

ist, welche die Abh&ngigkeit der Anzahl (N) 
der Punkte in einer Umgebung des betreffen- 
den Punktes von einer Abmessung (d) der Um- 55 
gebungdarsteQt; 

d) einer Vorrichtung (112) zum Bereitstellen 
von Skalierungskoeffizienten der Punkte einer 
n-dimensionaien stochastischen Punktevertei- 
lung eines der zu filternden Punkteverteilung 60 
aquivaJenten Typs; 

e) eine Vorrichtung (116) zum Selektieren der 
Punkte der zu filternden Punkteverteilung. die 
in einem Bereich der Skafierungskoeffizienten 
liegen. in denen sich die Skalierungskoeffizien- 65 
ten der zu filternden Punkteverteilung wesent- 
licfa von den Skalierungskoeffizienten der sto- 
chastischen Punkteverteilung unterscheiden, 
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und 

f) eine Vorrichtung zur Darstellung der Punk- 
te, denen die selekrierten SkaHerungskoefS- 
zienten zugeordnet and. 
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